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I. Introduc t ion 

The mechanisms of development of severe thunderstorms are not well 
understood. Rapid developments of storms and changes in atmospheric 
conditions can occur over small time and space scales, greatly influencing 
local weather and sometimes producing localized severe weather events. 

By providing atmospheric sounding data at 3-hour intervals rather than 
the usual twelve-hour Intervals, the AVE experiments allow finer time 
scale resolution of weather features and the possibility of more adequate 
understanding and prediction of processes. 

The objectives of this report are: 

1) To present a time series of maps of atmospheric parameters 
as measured for the AVE IV experiment. 

2) To present a thunderstorm model which incorporates mesoscale 
vortices. 

3) To explain a severe weather predictive index developed by 
Eagleman on the basis of a vortex model. 

4) To present results of the application of this index to the 
AVE IV data. 

5) To compare these results with the results from applying 
presently used predictive indices. 

6) To suggest modifications of the predictive index. 
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XX, Synoptic Conditions 

To determine the synoptic patterns of the thermodynamic and kinetic 
variables, surface maps and 500 mb charts were drawn for each time of 
the AVE IV rawinsonde flights of equivalent potential temppralure, water 
vapor mixing ratio, potential temperature^ and the magnitudes of the 
wind speed* 

Figu. 2 1 contains the surface maps and 500 mb charts for 0000 GMT, 
April 24, 1975. The surface map shows that much of the Central United 
States was covered with potentially warm, moist air. Values of water 
vapor mixing ratio in excess of 14 g 'kg were found in much of 
southeastern Texas and southern Louisiana. Mixing ratios of 10 and 
above were found throughout most of the Central Plains. Potential 
temperatures ranged between 295®K and 310®K across most of the experiment 
area. Maximum surface wind speeds were located in southern Texas 
and in a band across Oklahoma, central Missouri, southern Illinois, and 
Indiana, into West Virginia and Pennsylvania. 

The 500 mb chart shows the maximum values of mixing ratio extending 
through the eastern U.S. from Louisiana to Pennsylvania, Maximum 
values of potential temperature are also in this area. Maximum wind 
speeds are in West Virginia, Pennsylvania, and Ohio. 

Figure 2 showsvariable fields at 600 GMT, April 24, 1975. At the 
surface, large gradients of mixing ratio and equivalent potential 
temperature exist over Central Texas and Oklahoma and Kansas, with values 
of w varying from 4 to 16 across Texas. Potential temperatures over the 
area are 5^ to lO^K lower than at the previous sounding time. 

Wind speeds are also lower, with a local maximum at Topeka. 


2 



At the upper level, the area of maximum w has moved to the east 
and south. Potential temperature lines have moved southward, and the 
maximum wind speeds are over the Mississippi River basin. 

Figure 3 contains the surface map and 500 mb chart for 1200 GMT 
April 24, 1975. The large surface gradients of w and 9^ remain over 
western Texas, with an intrusion of dry, potentially cool air into 
northern Texas. Lines of potential temperature have moved further south 
and west. At 500 mb this cooling is also present. The location of 
maximum winds has changed little. The analysis of the area around 
Monette, Missouri, may not be a measure of synoptic conditions 
because the sonde was released during a rainstorm. 

Figure 4 shows the atmospheric conditions at 1500 GMT April 24. At 
the surface the intrusion of warm, dry air is still over northern Texas. 
Lines of potential temperature have moved northward in the eastern 
section of the country. Increased wind speeds are also evident. There 
is no feature at the 500 mb level corresponding to the high surface 
gradients of 0^ and w over Texas. Lines of 9 have moved north in the 
East and South. 

At the surface the potentially warm air has moved further into 
Texas. At 1800,. as shown in Figure 5, areas of high moisture have 
moved slightly eastward. 

At the 500 mb levels patterns are consistent with those of the last 
flight. The area of maximum 0 ^ has moved eastward from Indiana to 
West Virginia. Figure 6 contains the surface maps and 500 mb charts for 
2100 GMT April 24, 1975. A localized maximum of 0^ and w is present at 
the surface at Shreveport, Louisiana. Potentially warm air continues to 


move eastward across Texas, At the 500 mb level the moisture maximum has 

spread and is now centered over West Virginia and Kentucky. Figure 7 

contains the surface map and 500 mb chart for 0000 GMT April 25, 1975. 

The localize 1 maxima are no longer present in Louisiana at the surface. 

At the 500 mb level the maximum of 0 has moved over Tennessee, Kentucky, 

e 

and surrounding areas. Figure 8 shows the surface map and 500 mb chart for 
600 GMT, April 25, 1975. Strong gradients of equivalent potential temperature 
and mixing ratio occur in Texas at the surface. Potential temperatures 
are 5® to 10® K cooler over most of the experiment area than at the previous 
sounding. At the 500 mb level an area of relatively moist air is centered 
just west of the Mississippi River, and potentially cooler air has moved 
into the northern states. Figure 9 shows the surface map and 500 mb chart 
for 1200 GMT / pril 25, 1975. The patterns have changed little at the 
surface from the previous sounding. Potentially cooler air has moved into 
Texas. At the 500 mb level the area of maximum mixing ratio has moved 
eastward and the high velocity winds are centered over Middle Tennessee 
and Northern Alabama. 

A number of severe weather events of various types occurred during 
the AVE experiment. These include tornadoes, damaging winds, hailstorms, 
flash floodings, and funnel clouds. A complete description of the severe 
events is given by Turner [1976]. For the purpose of this report a 
severe event is defined to be a tornado that touches down, a wind with 
speeds over 50 miles per hour, or a hailstorm with stones greater than 
0.5 inches in diameter. A map of the locations and types of severe events 
is given in Figure 10. The reported severe hail events are shown on 
the north and west edges of the severe weather region and tornado and high 
winds are shown in the center, south and eastern portion. 
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Figure lb. Surface map of potential temperature and wind speed at 
0000 GMT April 24, 1975. 
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Figure Ic, 500 mb chart of equivalent potential temperature and 
water vapor mixing ratio at 0000 GMT April 24, 1975. 




Figure Id. 500 mb chart of potential temperature and wind speed 
at 0000 GMT April 24, 1975. 




Figure 2a. Surface map of equivalent potential temperature and water 
vapor mixing ratio at 600 GMT April 24, 1975. 
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Figure 2c* 500 mb chart of equivalent poten ial temperature and water 

vapor mixing ratio at 600 GMT April 24 , 1975. 






Figure 2d. 500 mb chart of potential temperature and wind speed at 

600 GMT April 24, 1975. 



















Figure 4a. Surface map of equivalent potential temperature cind water 
vapor mixing ratio at 1500 GMT April 24, 1975. 




Figure 4b, Surface map of potential temperature and wind speed at 
1500 GMT April 24 , 1975. 
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Figure 4d. 500 mb chart of potential temperature and wind speed at 

1500 GMT April 24, 1975. 




21 


Figure 5a. Surface map of equivalent potential temperature and water 
vapor mixing ratio at 1800 GMT April 24, 1975. 

















Figure 6c. 500 mb chart of equivalent potential temperature and 

water vapor mixing ratio at 2100 GMT April 24, 1975. 






Figure 6d. 500 mb chart of potential temperature and wind 

at 2100 GMT April 24, 1975. 







Figure 7b. Surface map of potential temperature and wind speed at 
0000 GMT April 25 , 1975. 
























Figure 9a. Surface map of equivalent potential temperature and water 
vapor mixing ratio at 1200 GMT April 25, 1975. 




Figure 9b. Surface map of potential temperature and wind sneed at 
1200 GMT April 25, 1975. 








Figure 9d. 500 mb chart of potential temperature and wind speed at 

1200 GMT April 25, 1975. 





Figure 10. Location and nature of severe events during the AVE IV 

experiment. H = hail greater than 0.5 inch in diameter; 
T = tornadoes; W = wind speeds in excess of 50 mph. 
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III* T he Tliuiiduf s Lorui Mode] 

The airflow in a thunderstorm is important to the mechanics of 
thunderstorm development. Connell (1973) has suggested that some 
thunderstorms might contain a pair of contrarotating leeside vortices 

based on aircraft measured winds near cloud base. The existence of these 
vortices has been substantiated by the dual-doppler radar measured winds 
obtained by Kropfli and Miller as shown in Figure 11. It has been 
suggested that a thunderstorm might act to block the environmental winds 
in much the same way that a solid cylinder would. However, it is known 
that a thunderstorm does not completely apprDximate a solid cylinder 
but entrains some environmental air. Using the experimental results 
of studies of jets in crossflows by Jordinson (1956) and by McAllister (1968), 
and the dual-doppler radar measurements of Kropfli and Miller (1975) , Connell 
(1976) has proposed that the interaction between a thunderstorm and its 
environment might be analogous to the behavior of a jet in a crossflow. 

In order to produce a blocking effect in the cloud layer environmental 
winds, the inflowing updraft must be opposed in direction to some of the air flow 
above the cloud base. Since the thunderstorm itself usually moves, it is 
the relative updraft direction that must oppose the relative environmental 
winds . 

Thunderstorms usually develop in environments with strong wind shear; 
and, in fact, the double vortex model requires some strongly sheared 
environment. As the low-level updraft enters the clouds, it encounters 
upper-level winds from an opposing direction. The intrusion of the 
updraft into the upper level wind will produce a blocking effect and 
environmental air will flow around it with some momentum mixing which 
causes the updraft to bend downwind. Flow around the updraft core will 
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produce a lee-side wake where contrarotating vortices develop. Figure 12 
is a schematic of a double vortex thunderstorm. The intensity and 
vertical extent of these vortices is in part dependent upon the extent 
to which the cloud layer winds oppose the updraft inflow. 

A review of the literature ,has revealed a thunderstorm model 
developed by R- Eagleman (Eagleman, 1975) that is very similar to the one 
just described. Eagleman had used his model to develop an index for 
predicting the occurrences of tornadoes. This predictive scheme, 
which will be described in the next chapter, was applied to the AVE 
IV data^ and results are given in Chapters V and VI. 
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Y Thunderstorm 



a. Horizontal flow relative to 
moving storm at 6.4 km above 
ground. Lengths at arrows are 
proportional to relative wind speed. 


Thunderstorm 


X=5.4 



0 6 12 18 
Y km 


b. Vertical flow in y,z plane 
at X “ 5.4 km. 



c. Vertical flow in x,z plane 
at y • 11.4 as indicated in Figure 
11a. 


^ Storm Motion 



d. Schematic representation of 
flow pattern in vertical plane 
corresponding to Figure 11c. 


Figure 11. Three-dimensional velocity structure within a severe 
thunderstorm measured by dual-doppler radar (from 
Kropfli and Miller) . 





IV. The Predictive Indices 


The accurate forecasting of tornadoes is of importance to the 
protection of life and property. At the present time forecasting techniques 
are not accurate. Today a tornado will occur forty percent of the times 
the National Severe Storms Forecasting Center issues a tornado forecast. 
Furthermore, only thirty-five percent of tornadoes that do occur are 
within the tornado watch area. The need for more accurate predictions 
is evident. 

A predictive index has been devised by J.R. Eagleman which incorporates 
both atmospheric thermodynamics and environmental wind conditions hypothetically 
leading to double -vortex thunderstorm formation. Darkow’s Energy Index is 
used as an indicator of the potential instability of the atmosphere. It 
is combined with a Shear Index, which reflects the shear in the relative 
winds between the surf ace-to- 850 mb layer and the cloud layer, to produce 
an Energy Shear Index. Each of these indices will be explained in the 
remainder of this chapter. The Severe Weather Threat Index (SWEAT) and 
the Surface Potential Index (SPOT) , which were also calculated for the 
AVE IV data for comparison with the Energy Shear Index, will also be 
explained in this chapter. 

1) The Energy Index (El) 

The total specific energy of a mass of air may be expressed as 
2 

E(T)= CpT + gZ + Lq + V /2, the sum of specific enthaply, potential 
energy, latent energy, and kinetic energy, where c^ is the specific heat 
of air at constant pressure, T the temperature, gZ the geopotential, L the 
latent heat, q the specific humidity, and V the scalar velocity of the wind. 




Darkow (1968) points out that the kinetic energy term is usually two 

orders of magnitude smaller than the other terms and may be neglected. 

The errors in reported upper air humidity values allow 

the additional approximations q = w and L = where w is the mixing 

ratio and a constant latent heat of condensation. The total energy 

or static energy Is then expressed as E(T) = C^T = gZ + total 

energy defined by this equation is conservative with respect to both 

unsaturated and saturated adiabatic processes. Thus, the potential 

convective instability of atmospheric layers is indicated by the amount 

of decrease of total energy in the layers. Using the values of 

c a 1.00 J -2500 J and g = 980 cm sec"^ yields E(T)= 

p g ’ o 8 

-3 -1 

T(K) + 9.8 x 10 Z(m) +2.5 w(g k^ ) which may be approximated as 
E(T) =^T(®K) + 2.5 w (gkg'^) + Z(m)/100. 

A stability index, the Energy Index, whic!>. reflects the contribution 
to total atmospheric energy of both ascending potentially warm air and 
descending potentially cold air Is defined as the algebraic difference 
between the total energy of the air at the 500 and 850 mb levels, 
expressed as El = E(T)t.rtn“ E(T)„cr>* This difference is shown schematically 

jOU ODU 

in Figure 13. The 850 and 500 mb values are chosen as representative of the 
low-level air and mid-tropospheric air entering the storm as dictated 
by the routine availability of data at these levels. 

Quantitative values of the Energy Index had been assigned to various 
degrees of thunderstorm intensity. In the range of 0.0 to -1.0, 
thunderstorms are possible but will not be severe. In the range of -1.0 
to -2.0, isolated severe thunderstorm activity is possible, particularly 
as a continuation of severe activity moving into the regions. For 
Energy Index values less than -2.0, severe thunderstorms and tornado 
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activity are highly probable, providing an adequate triggering mechanism 
to release the potential instability is present. 

Shear Index 

One triggering mechanism that is present in the environment of a 
cloud is the vertical shear in the relative winds between the low levels 
and the mid-troposphere. In this regard, Eagleman (1975) has developed 
a Shear Index to reflect vertical changes in relative vector velocity. 
Figure 14 shows the relationships between cloud motion, measured winds 
relative to earth, and the winds measured relative to the 
moving cloud. The mean environmental wind is calculated by finding 
the vector mean of the 850, 700, 500, and 300 mb reported winds. These 
are interpolated sc that the relative velocity through each 50 mb layer 
can be calculated. 

In order to calculate relative velocities, storm speeds and 
directions must be known. For the calculation of the Shear Index, seven 
storm speeds of 50, 55, 60, 65, 70, 75 and 80 percent of the mean 
environment wind speed and 26 storm directions ranging from 60 degrees 
to the left to 60 degrees to the right of the mean environmental wind 
direction, incremented every 5®, were used. This produced 182 variations 
of storm movement as shown in Figure 15. 

Since the model requires blocking of the inflowing low level air 
in the mid-level, layers 150 mb deep are examined to see if they oppose 
the relative winds of the low level. 
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The surface-8'jO nit layer Is usi‘d as the low level inflow layer. Six 
layers were chosen as the critical mid-level layers for the occurence of 
opposing component velocities^ these were the 650-500 mb, the 400-250 mb, 
the 550-400 mb, the 500-350 mb, the 450-300 mb, and the 400-250 mb layers. 

The opposing components in the six mid-level layers were considered to 
produce blocking if they were between 75 and 125 percent of the magnitude 
of the surface-850 mb wind. 

Thus the Shear Index (SI) for a trial storm and speed was defined 
as the number of consecutive mid-level layers whose opposing components 
were between 75 and 125 percent of the surface-850 mb wind speed. Therefore, 
the Shear Index can vary from zero to six; a maximum Shear Index of 6 
is shown in Figure 16. The Shear Index of a sounding was defined as the 
largest SI for all 182 trials obtained for the 7 storm speeds and 26 
directions. The direction of movement of an actual storm should correspond 
to one of the calculated storm directions which yields the maximum SI. 

The Energy Shear Index 

The Shear Index measures only the atmospheric wind profile; the 
proper thermodynamics must be present also for storm development. 

Therefore, the Shear Index is combined with the Energy Index to produce 
the Energy-Shear Index. To determine the best empirical combination of 
the Indices, the SI was graphed versus the El for 59 soundings as shown 
in Figure 17. Twenty-seven of these soundings were proximity soundings; 
twelve were precedence soundings, and twenty were nonproximity soundings. 

A proximity sounding was defined as a sounding within the warm air sector 
and less than 120 miles from a coniirmed tornado touchdown, and within 
two hours before the tornado or no more than one-half hour after the 
first report of a tornado. A nonproximity sounding is for the same 
time period but located over two hundred miles away from a tornado 
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occurrence. Precedence soundings are tho^e taker in the warm air ahead 
of the cold front but removed from it in either time or distance. 

A line can be drawn that separates the proximity and precedence 
sounding from the nonproximity soundings. The equation of this line is 
El = 1/2 SI-2 or 2EI -SI + 4= 0. The Energy Shear Index can thus be 
calculated from the equation ESI = 4-SI + 2EI. If the ESI is negative 
and a cold front id nearby, tornadoes are predicted; if it is positive, 
tornadoes are not indicated. For a more detailed description of the 
Shear Index and the Energy-Shear Index, see Eagleman (1975). 

Sweat and Spot Indices 

As a basis of a comparison of the accuracy of the Energy-Shear 

Index, two currently used indices, the Severe Weather Threat Index 

(SWEAT) and the Surface Potential Index (SPOT) were calculated. These 

indices are used in conjunction to produce short-term (three to six hour) 

depictions of areas with high potential for severe storm development or 

occurrence. The soft SWEAT was calculated in this study since parameters 

from AFGWe Fine Mesh and Boundary Layer Models were not available. 

The Sweat Index is computed using the equation SWEAT = 12 t^+ 20 (T-49) 

+ 2W + W^^^ + 125 f/2) where 
e juu 

t - low level dew point in "C, the level used is 850 mb 
e 

in the soft SWEAT and 900 meters in the BLM 
computations. 

T = Total totals (T=850 mb temperature plus 850 

mb dew point temperature minus twice the 500 mb 
temperature, all °C);for complete details on this 
stability index see Miller (1972). 
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W^*= low level wind speed in knots; the level used is once 

again either 850 mb or 900 meters. 

500 mb wind speed in knots. 

f(a) = a step function of the veering angle to 

See Figure 18a for a plot of this function. 

This term is set to zero if both W and are not 

e 500 

equal to o:* greater than 15 knots. The term is not 
computed unless the 850 mb wind direction is within the 
range 130° to 250° and the 500 mb wind direction is 
within the range 210° to 310°. 

All negative terms are set to zero. 

The SPOT index is computed from the equation SPOl = (t-60) + (t ,-55) 

a 

+ 100(30.00 - p) + f(v) where 

t = surface temperature in °F. 

t,= surface dew point in ®F. 
a 

P = altimeter setting in inches. 

f(v)= wind speed term which is determined from tht» table shoKn; 
in Figure 18b. 

Negative values are allowed to occur. The altimeter term is reduced by 
50 percent when temperatures are less than 50°F and altimeter settings are 
below 29.50 inches. Regions where high values of the index lie in close 
proximity to very low values are suspect areas. 
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Figure 13. Schematic total energy profile in late afternoon prior 
to the outbreak of severe thunderstorm activity (from 
Darkow, 1968). 





Figure 14. Vectors for a tornado proximity sounding showing the 

wind (R) relative to a moving thunderstorm as determined 
by the movement of the storm (t) , which creates a wind (i'^) 
opposite to the direction of movement of the storm. The 
combination of (r ) with the measured winds (m) relative 
to a fixed point results in the relative winds (R) for a 
moving thunderstorm (from Eagleman, 1975) . 
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Figure 15. Combinations of storm speeds and storm directions used for 
making calculations of relative winds and shear index 
(from Eagleman, 1975) . 
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Figure 17, Scatter diagram of SI versus El used in determination of 
equation for ESI (from Eagleman, 1975 ) . 
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Figure 18a. Step function used to determine f(2) term in SWEAT 
equation (from Miller and Maddox, 1975) . 
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Figure 18b. Table used to determine f(v) term in SPOT equation 
(from Miller and Maddox, 1975) . 
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V. Applications of Indices to AVE IV Data 

A computer program was obtained from Dr, Eagleman for computing the 
Energy Index, the Shear Index, and tlie Energy-Sh'»:ir Index. The input 
format parameters were changed to read the AVE IV data, and the section 
of the program which determines the atmospheric variables at 50 mb intervals 
by linear interpolation was removed since data at these intervals are 

directly available from the AVE IV data. The program as received did 

not execute properly on the IBM 360. For certain angles between the 

assumed storm direction and the measured winds, division by zero was 

produced. Also the accumulation of computer-generated round-off errors 
sometimes produced values of the sine or cosine of an angle whose 
absolute value was greater than one. After these problems were eliminated, 
a subr*)utine was added to calculate the SWEAT and SPOT indices, and the 
program was run with the AVE IV data from 29 stations at all nine times. 
Figures 19 through 54 show the results of these computations. The 
shaded areas on the maps of Energy Index are those areas in which El 
has a value more negative than -2 and severe thunderstorms and tornadoes 
are probable. The weather types associated with other values of El have 
been explained in Chapter IV. Shaded areas on the ESI maps indicate those 
regions with ESI values less than zero. Severe event locations and types 
are indicated. The time of the sounding was between two hours before 
and one hour after any severe events marked at a sounding time. 

The following comments summarize the tests using AVE IV data: 
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A summary of the number of occurrences of each type of severe weather 
is given in Table 1, along with the number that was correctly 
predicted by each index. The sweat and spot index predictions are not 
included because it was felt that the criteria for both are not 
sufficiently definitive to permit an objective prediction. As we gain 
experience with their use we will better trust our subjective use of 
these indices and may use them for comparison. 
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Table 1 


Events Reported 


Events Correctly Predicted by 

Type _ . _ 

Number 

El 

ESI 

Hail 

4 

4 

4 

Tornadoes 

9 

4 

6 

Wind 

5 

4 

3 

TOTAL 

18 

12 

13 
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Figure 19. Map at 0000 GMT April 24, 1975 for Energy Index. 








Figure 22. Map at 0000 GMT April 24, 1975 for Sweat and Spot indices. 





Erergy Index. 






0600 GMT April 24, ’975 for Shear Index. 








Figure 25. Map at 0600 GMT April 24, 1975 for Energy Shear Index, 







Figure 27. Map at 1200 GMT April 24, 1975 for Energy Inde 




Figure 28* Map at 1200 GMT April 24, 1975 for Shear Index. 
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Figure 38, Map at 1800 GMT April 24, 1975 for ^weat and Spot indices. 





Figure 39. Map at 2100 GMT April 24, 1975 for Energy Index. 








Figure 42. Map at 2100 GMT April 24, 1975 for Sweat and Spot indices. 









Figure 45. Map for 0000 GMT April 25, 1975 for Energy Shear Index. 
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Figure 48. Map at 0600 CMT April 25, 1975 for Shear Index. 








Figure 51. Map at 1200 GMT April 25, 1975 for Energy ":*Jex. 



Figure 52. Map at 1200 GMT April 25, 1975 for Shear Index. 







Figure 53. Map at 1200 GMT April 25, 1975 for Energy Shear Index. 









VI. 


Summary of Results 


A graph of the El's computed in this study versus the 

Si's is shown in Figure 55. While nine of the twelve 

proximity soundings and all but one of the precedence soundings fall below 
the line El =0.5 SI-2, this line does not uniquely separate these 
soundings from the nonproximity soundings as it did for the soundings 
studied by Eagleman. (See Figure 17, Chapter IV, for comparison,) A linear 
relation between an ESI and si and El is not clearly indicated. 

However, the tests of the ESI . computed on this basis do show that the 

ESI so computed does improve the forecast of severe wind events due to 

thunderstorms (such as tornado and linear winds). It is slightly better 

than just the energy index which has no windshear component. It also 

is better than the SWEAT which does Include an average wind shear related factor. 

The results with the small sample from AVE IV are not as striking as 

those shown by Eagleman. (See Chapter IV Reference, Eagleman, et al. 1975.) 

VII . Conclusions and Recommendations 

The incorporation of vertical shear of horizontal environmental wind 
(1) in smaller layers ("200 mb thick) and (2) in a manner consistent with 
a double vortex model of a thunderstorm does produce slightly Improved 
prediction of /or positive correlation with severity of thunderstorm winds 
and velocity of thunderstorm motion. It is still true that many 
thunderstorm severe events and velocities of movement are not predicted 
by this or any current method and the false alarm rate is high. 

The method of Eagleman requires making a large number of guesses 
of the direction and magnitude of the motion of the predicted storm. 

Further, it does not prevent prediction of several directions of motion 
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at the same time. Finally, the choices of positions of significance 
for voftlces and the thickness of layers for computing wind shear are 
neither as refined nor as physically based as now seems desirable. 

It Is recommended that steps be taken to Improve physical Insight 
Into the mechanisms set Into operation by the Interaction between the 
thunderstorm and the sheared environment. However, even without the 
detailed Understanding of these processes It Is possible to test an 
Improved Shear Index. Several suggestions follow: 

(1) Use thinner layers of atmosphere for computing shears. 

(2) Reduce the amount of guessing of storm vector velocity which 
Is required in the present method. 

(3) Improve the form of the energy shear index relation. 

(4) Eliminate cloud zones of probable irrelevance from the shear 
index calculation scheme. 

(5) Utilize more field experiment data both for developing better 
correlations among prediction parameters and for testing the 
predictive schemes. 
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Appendix A 




!• A printout of the computer printout is attached as exhibit A-i. 

2. Program Checkout. 

To insure that the Energy-Shear program was making the desired 
calculations for the shear index correctly, the program was run with 
the data from Charleston, South Carolina, at 2100 GMT on April 24, 1975. 

An assumed storm speed and direction of 50% of the cloud layer average 
winds peed and 10 degrees to the left of the average wind were chosen 
for a test calculation. These parameters gave a computed storm velocity 
of 5.45 meters sec ^ from 274.63 degrees. The shear index was 5 with 
surface to 850 mb average winds of 8.10 meter sec ^ from 179.97 degrees 
and a 550 to 300 mb average wind of 9.92 meters sec ^ from 315.47 
degrees. 

Instead of following the computer program step by step, vectors 
were used to calculate the relative mean wind velocities for each 150 mb 
layers and the component of these velocities opposed to the surface-850 mb 
mean velocity. 

The components of the cloud layer average wind were calculated by 

finding the average of the u and v wind components at the 850,700, 500, 

and 300 mb levels as given in Fucik and Turner (1975 ). The average winds 

for each 150 mb layer were also found using the published u and v values. 

Table 1 gives the components of the mean wind and the average wind 

through each layer, and the 550-300 mb average wind. 

- -2 -2 ^ 

The mean cloud layer wind was calculated from ^ ^cl^ 

the direction from 6« Arctan (=^ — ). This yielded V -■ 10.8783 and 

u ^ cl 

cl 

6 <■ -14.643. It must be remembered that this 6 is the angle in 


the usual mathematical coordinate system and must be converted to the 
meteorological system. 

Remembering that the assumed storm speed and direction correspond 

to 50% of the cloud layer average speed and move 10 degrees to the 

• -1 

left, the assumed storm moves at 5.44 meters sec from 274.64 degrees. 

This velocity was also resolved into u and v components to facilitate 

calculation of the relative mean wind in each layer. Relative velocities 

were calculated using the equation ■ Vj^. The components of the 

relative wind are given by Table 2. 

The component of the relative wind in the upper layers which directly 

opposed the surface-850 mb relative wind was then found by 

V „= • “ - , r ~ — . This component is also given in Table 2. 

U-B Vg 

The magnitude and direction of the surface to 850 mb layer relative 
mean wind were calculated from the components by the method described 
earlier for the cloud layer average wind. The resultant relative wind 
was 8.1050 meter sec~^ from 182.04 degrees. 

Since the model assumes that the component of the relative mean wind 
velocity of an upper layer that directly opposes the surface-850 mb 
relative wind must be between 75 and 125% of the magnitude of that wind 
in order to produce blocking, the value of V^_g should be between 
6.075 and 10.125 if the layer is to have a blocking effect. From Table 2 
five consecutive layers have values of within this range; therefore 

the shear index is 5. Table 3 gives a summary of the values computed from 
the IBM 360 and those calculated using vectors and the HP-65 calculator. 
Discrepancies between these values are probably due to sounding differences 
in the two machines and are well within the accuracies of the data. 
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Exhibit A-1. The Computer Program 




- s JCB 


T1MC«^#PAGCS«99 


C TM|G-P«OG«AM calculates ThE E«i£*tGY«Sf»€AP- INDEX F«»ON 9AwfN5ON0t -DATA 

C 

c 

— c .... 

1 DI^CNS10.N 3Trvi3)»0ATEI9)# P^ESSi 3) #6P-iT i 3) •Tev»( 31 *£ T C 2 i t UU 

2 3r<E;SlON WDUDf DIR MDC6t2'J0l « (6 t2U0) tRuS ( ) 

. 3 .-- . 01«EN3|CN CwS«A0)tJ/5523«A*2O'>)f •*SD29*. 6 •2C3I t DEV23f6*2CD) 

4 dimension SPAS23C6>2CC}fSRAD23(6«2:.MSS23t&t£:.d)tSD23(9*20J) 

5 DIMENSION ISAVE(6» 

6 INTEGER HOUR * .... 

C 

7 DATA AuEFT/4HL€FT/ 

8 . - DATA RIGNT/4MRGHT/. 

C 

9 COMMON TE^'P* OrCMPt wOt 4St HOUR* NOAY* PRESS# GPhT 

13 COMM ON I ST A 

C READ 1\ A DATA SET 
C 

.11 .. 100 READ CS#1102) ISTa# DATE - - . 

12 IF I ISTA .EC. 999991 STOP 

13 read (SfllOL) GPHTm* »'R£S5 U»f TEmP U)» DTEvod!# UDil >• ^Sl 
1 HtH0JR»\3AY 

14 03 lOS L0P1» 2f7 

15 READ (SfllOS) w? (LOPl)f WSILOPII 

U... lOS CONTINUE - . > - . - - 

17 READ :S#11C4I CPHTdIt PRESS (2)# TEVP (2lt 0T£/P(2)» wDia )• 4St 
I 8) 

18 S3 110 L0P2* 9»21 _ . ~ - 

13 READ './3 (L0P2I# N$(L0P2' 

20 110 continue 

2i. . PEAO CS#1104l 6PHTI3U PRESS I3U 7 CmP (3»« or£MP43)» wO(22)t WS(- 

1 221 

22 03 IIS LCP3* 23t41 

23 READ (S#113S1 WO (LOPIlt WSILOPSl ... 

24 119 continue 

25 no? format (1S»2FX«SA4I 

28 _ 1104 format (2r7df1PA*l«FS.l»3SXflS«13t - — 

27 IlOA format (2AX#F6«lt^S#l» 

C IF DATA IS MISSING# PRINT veisAGE 

28 IF t SPhT(2) «>’3# 09#9 .OR# OPhTISI .ED. 99.9> GO TO 299 

29 1* I TEMP(2I .EQ. 99.9 .OR. T£moi3) *£0, 99«9) GO TO 299 

' 30 - IF t DTEMPI2» #E0# 99#9 #0R#DTFMPI3! .EO# 99«9) CD TO 299 


C 



c 

C calculate The energy index 

c 

31 oe 240 I >2 #3 

33 tCM-'A « DTF.4P lit ♦|73#l8 

31 TEv^tt • TCv?i in ♦ZTS.U 

34 |F |T«MPA .LT. GO TO 22^ 

35 EXP • II7,2693M2 «( TE 4VA*273#I8M / tT|MPA*|!#«88» 

36 GO TO 233 

37 ??C EXP • »’UX^4118*| Tfv.»A»273.l8nnTEMl»A*7.88l 

38 213 r • 8tn7At(7.TlA.^A**lXPt 

10 A*'tx;9A • t622#«EI/PH:ESSl It*Et 

40 J*l«l 


> 

1 

u 

♦* 

r 




is 
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ETU» • .24»{TEMPft ♦2.54AVIX^A4CPHTm/X00.l 
240 CvONTiNUr 

El- ET(?)-ET(1) 

CO TO 3C0 — 

299 W9ITE (6«1012) OATEt ISTA 

10X2 F09PAT (42HTHE ENERGY INDEX CANNOT BE CAtCULATED FOR tBA4»XX« 
- GO TO 100 


: calculate the mean wind velocity of tme clcuo layer 
300 4RITE (6«XX22) DATE# IfTA 

1122 F0RMATtXHX*///»X7X#SA4flXfI5*40*»X8H$HEAR INDEX VI MOO#//) 
tCI- 0 

IFOOR • 0 - ~ - — 

DO ?04 I-l#20 

IM.i 

IF n .NE* X) IM-2AI-2 - 

IF (itfSMM)«NE» 99-9) GO TO 304 
IC1> ICI4X 

. IF ClM.rc.d «OR. IM«C0»14 »09t I*^*EQ»22 *OR# I^#EO«39> GO 

GO TO 304 

302 IFOUR • IFOUR41 

304 CONTINUE — 

304 N ■ 20«rci 
XWINO «C# 

. . Y;«IND « 0»- - — - 

00 3X2 j«lf20 

jM.j 

IF ij.Nr. n UM> 2 «J «2 .... - - - --- 

IF US(JM) •EOt 99«9I 00 TO 312 

IF »UR# ^•*<«EO«22 -OR- GO 

00 TO 112 - - 

310 3A0«^S ■ -0UM)/57«294 

XWIN3 ■ XKlND « 4$UM)«C0SIRAMSI 

t»{lND ■ th;l\D ♦ *$IUFtl*SlN4RA0«SI 

3X2 CONTlfiUF 

320 AVXW « XWlNO/|4«-FLOATUFOURn 
AVYW • Y«|ND/(4**FL0AT(IF0URII 
TM<S« 509TIAVXW4»2 4AVVW442I 
DO ■ AVYW/AVXW 

TMa'D • AT4N(D0)*3T«294 * — - 

|FIAVXx,C£#0») 00 TO 390 
Ty#»D ■ TMaO ♦l§0* 

3%0 IF tTvwD *LT, 0«) TvwO • Tm<#D4360# 

IF ( TVaD *3T* 340# ) TMWD ■ TMaO-360# 

DO 363 
ISAVfIVl «X 
363 continue 
INDSH3 • 0 
DO 699 Ills 1«3 
DO 699 X - l#6l«3 
DO 699 L • l#31f« 


C 

C calculate relative aind velocities BASFD on deviations 

C OF STOR^: DINECTIONS FNOM ThE ^EAN «|\D VELOCITY 
C 


■uINAL 

OF wu 





to IMUNWI* 50^L-1 

91 ICEVIT • <-l 

9? SS • Tvw$ •FLOAT (90^fl)/100* 

9Y — .. If |llt«£0«l> 019* fllOHT ♦ • 

94» IF (HUES* 2\ 0IR« ALEFT 

99 ir I DIR«FQ» ALFFTI lOEVlT* lOEvlUf 

96 - -- IF I II1«E0«2> GO TO 401 - - ^ 

97 SO* TMwO ♦FLOAT UJ-W 

90 GO TO 402 

99 401 SO* Tm^<d»FLOAT{ FI-4, — - -- ♦- 

IOC 402 DO 490 LL*1«2C 
XOl LLV • LL 

102 -- IF <LL«SE«1 lLLM*2^Ll.-2 - -- - — 

109 GAvva •SO-»iD(H.VJ 

104 9A0GA* GAvyA/97*296 

1®9 ' AwSlLLfc*! • S09T(WS(LLM)#*2^SS«^2«2*0«4S(Lt.vt«SS*C3S(9A9GAf ) * - • 

106 SI\SS * (SS*SIN(RADGA| )/R4;SILLMI 

107 N5IN «3 

100 IF (SINSS,LT,0*0I NSIN*1 - . . . 

109 IF (A«Sf5INSS».G£. .9999990 .AND, ABS(S!\SS » .lT . l.OGCCl ) Si\SS*l.O 

UO IF tS?\5S.GT.l. » 30 TO 1224 

111— • IF (\$I^ ,E0,1 .ANOt SINSS *F0« 1,0) SlNSS • -1«0 
112 IF ni,C-A2S(SlNSS)**2.I •GT, 0.0) GO TO 409 

119 THIS, 99021 III tK,L.LLM,Sl\S$,COSSS 

114 coss$*o,o 

119 CO TO 410 

116 409 COSSS •SQRTll.O- A0SI StNSS)*»2«3) 

117 9902 FOQVAT i4l4*2EU.6l - 

llA 410 'F (COSSS .LEt ,0000011 GO TO 4l| 

119 TANSS* sinss/cosss 

120 nCTA « ATAMTANSS)*97,296 

!?: uO TO 42C 

122 419 49ITE l6f99C2) 1 1 1 tKvLtlLMfSX^SStCOSSS 

129— IF ISINSS - 0,0) %16,4U,417 • - 

124 416 SETA • 270* 

129 GO TO 420 

126 617 BETA * 90* 

127 429 9W0ILLM) *uD (LLM)-9(TA 

120 IF {RwDlLL^ULr,0,l RwOCLL^I* 2*DaLKUI60, 

129 . IF IR40ILLMUGT«160* I AtfOUiMI* B^OlLLMt-lAO* 

190 490 CO\riNUE 

191 191 • 0 

192 DO 440 V* 1*S 

199 MMtV 

194 I* i» MM«2*y-2 

199^ IF 99.9i GO TO 440 

196 ici * n:^i 

197 440 CONTINUE 

190 AS'JV * C, 

199 VSUM * 0, 

140 NN * 9-I0I 

141 00 490 n*l«9 

147 nv.M 

149 IF I tI,‘4C,;i IIMv 24IX-2 

144 IF I wsdlvurc, 99,91 GO TO 400 

l49 990 • 9*r>| tlN*l/97,296 

146 ASjv, t |^'I*CDS M'4 xOI 

14? V$UV * vsuv#*,ftSI lly)4SI4IARADI 

149 490 CONTINUE 

c 


% 
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c 

C calculate I.AVEH AVERAGES 
C 


149 


AVXSJ'* • *SU^/6lOATINN> 

190 


AVV$JV ■ Y$uy/6L0ATCNM) 

^ 191 


AV^*S • SORT! AVA9uy«^2^AVVSgM442>- 

192 


0 • Avvsu'»/Avx$yy 

199 


AVOMrO ■ ATA*«C0)«97«296 

- 194 


|6 «AV*9U^«0^* Oil GO TO 901 ” • *' — " 

199 


AVRV.3 • AVO’WD4l60« 

196 

991 

IF ( Avt^*D iLT. Oil AV9/.0 • A^*«a0^960« 

197 


IF (AV^^O iGT# 360# > AyRwO • AVR40»960#0 ' 

199 


CVO ■ AVny^D ♦180# 

199 


IF f Ca uor# 360#} C40 • CWD-960# 

160 


00 970 «• 16»98f? - 

161 


<KJ «0 

162 


X9W *0# 

169 


• 0# - 

164 


KKl ■ »flC-6 

169 



- 166 


00 970 6<c« <FI#«K#2 

167 


(»Sl«Ll#rO,99#9) GO TO 910 

169 


999^0 •9x0I«L1/97.296 

169 


s X.^A^6«St<ALI*CCS(M240> - - - 

170 


v9m • ♦.^a$IAKU*S1\(RR9«0I 

m 


60 TO 920 

172 

910 

CiJ* KtJ^l - - 

179 

970 

continue 

174 

92* 

ir IKI(J«E0#4) GO TO 990 

179 


6AV9*' • A<i*/l4#«FLCAt UOII 

176 


V4V7# * Y7 f F<J) ♦ 

177 


aAV9«S ■ 60<4TUAv^a**A* TAV9w«»2> 


. 

e« VAV94/6AV^ii( - — * _ ^ 

179 


XAV9«K0 • ATAME)497.296 

190 


|F t XAV9i^#GE#0#) 00 TO 990 

161 


XAVIaO • AAV9«0«180# 

191 

970 

IF UAV<‘i»0#LT#0#l XA^9»0 • XAV9 f0«960« 

199 


IF 1 XAV940#GT# 960tl XAV9*D • XAV940«960# 

194 


CniGl ■ IC«0«XAV9WU> 

199 


IF (C0SG1#GT#90#I GO TO 960 

196 


CCEC2 ■ COECl « 90# 

197 


tr (COEG2#tT#a#l GO TO 960 

196 

940 

COCO • C0r:01/97.296 

199 


C«SUK) #XA/XF9tC0S ICOCG) 

— 160 

' — 

GO 70 970 

19: 

9*0 

C^SKO 

192 


GO 70 970 

>99 

960 

XAV9«0 • XAV9MD ♦ 960« 

194 


COCGl • XAV9W0 *CaD 

199 


XAV^tiO • XAV'IWV960« 

196 


16 ICr>fGl#L7. 47,1 GO TO 940 

197 


C*xD • M ♦IhO# 

• 199 


COfGl • C40 *XAV940 

199 


Q'MZ • C*C«963* 

200 


IF icorci #CT# 90# 1 GO 70 940 

201 


CWS U9l4 «!• 

292 

170 

CONTI Hilt 

^ 209 


TMIO 40# 

204 


XMIC 


..u'4AL PACr3 IS 

■b’ POOR O’’ 


20 % 

207 
20« 

209 

210 
211 - 
212 
21) 

214 - 
219 
216 
217 

219 

219 

220 - 
221 
222 
223 
226 
229 
226 

C 
C' 

c 

C CAuCULATf SHE69 |N0E< 
C 


227 

600 

IC29C 40. - - , . - 

229 


00 611 UK • 16f20t2 

229 


17 ICmSUJIU .LT. 9.) 60 TO 611 

290 


01F7 ■ C«6( 

2)1 


t::6» *CT. 0.1 33 TO 670 

2)2 

611 

CUNTINUt 

299 . 


DO 61» 1X0 • 22.92.2 — • 

214 


If (Ci^S (IFJ) .OT. 9.) 60 TO 620 

239 

619 

cosriNuc 

296 


GO TO 670 . .... 

297 

629 

6- UJ 

296 


6.1. 

299 


IXC29C .0 • • * 

240 


00 669 JJK* I4f92.2 

241 


17 IC4S IJJX) .LT.O.I 60 TO 629 

242 


Ot'F a CaS(JJX) •AV9k6 

249 


AAC7 a .296AV)X5 

244 


17 IA6$(0I7ri,LC. A6C2) 60 TO 699 

269 - 

679 

6 .'O. 

246 


30 TO 669 

747 

6)^ 

16 1 F .to. 0 .» 60 ro 663 

166 


IC29C • 1C29C ♦! 

269 


60 TO 690 

290 

66C 

IC29C •! 

291 

690 

16 1 1C29C .CT. 1XC29CI 1XC29C • IC29C 

2)2 


• • U 

29) 


30 TO 669 

294 

669 

CONTINUE 

219 


16 1 1XC79C .3T. t\OSH<) 1N0&H9 v IXClIC 

216 


|6 1 MC79C .10. 0» 60 TO 670 

297 


6C TO 6X0 

2)6 

673 

16 llll.ro. 1 .AN). X.CU.l .and. I.IU.II 60 TO 690 

|69 


CO TO 699 


C50 960 JJj«20f30»2 
ir CWSIJJJ) «E;0* 99,9) GO TO 979 

• 6WOUJJ//97.296 

xv|0 • K^'IO ♦9^SUJJ)#C0S CR9fl*9l 
Yv;0 • 7MI0 ♦H»S« 

60 TO 960 • - 

979 

990 CONTIM^E 

16 (JWJJ.JO* 6) 60 TO 990 — - 

XAVVJO • XVI0/(6,.fl0ATUJJJM 
V6VVI0 > rvro/(6*-6L5AT|JjJJU 
AVM«;PD • SwRTiXAV^lD ••2^ V4VM|0*92) 

C • VAVVI3/XXVM0 
Av'^OI^ • *TAN(0»*97.296 

16 lXAV«^IO#C.Ei'>*> 60 TO 999 — 

AVMOI<? • IVV019 ♦190. 

9»9 16 t AVM0I5 .LT.O.) AWDIR • AWOJX ♦360. 

16 ( AW60Z4 360.) AVMOIR • AvM>Z<4-36C. 

CO TC 6CC 
9«D Av**SPD • •!, 

AVM0I9 • *1. , 
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r 


6^0 • * 

a iswriixC2^c> 

vSSPi .UC25C*^V» • |mjnP« 

«50P5 <UC25Cf^^» • lOCVlT 
DFV2'. <UC25C»^M> • Din 
$RaS 2^ i|xC2iC*wM» ■ AVRWS 
Sfl*<02» <UC25C#v'^» • AVftifO 
S$2iaxC?9C#vv» • ss 
S025 UJ(C7SC**^^) -SO 
0|B''IO MxC75C*»»'»» • AV?^IH ■ 
5»:viOeUC25C#-^'') • AVV.SJ»0 
ISAVF( IXC25C) • I$AVCItxC2SC) ♦! 
GO 70 699 
69C I2E90 «0 

isxvEs • iMuspn 

ISXVEO • lOfVlT 
SAVCDV • DI9 
SAv£«($ • AV9iilS 

SAv£m<D • AvnwD 
SAVES *5$ 

SAVED* SO 

SAvE'^ • AV^Din - - . - 

f.AvEM$ ■. AVMSPO 
699 C^^T1^UE 


CALCUkAtr CNE(iGY-SMEA9 INDEX 


ESI * Aa -float (|NDSh9» ♦ 2« •£! 
IF U2E30.CS.0) GO TO 920- 
00 I ■W<‘ 

r I ISAvCiifaEOa U GO TO 9id 
HMlfm ISAVE 1 1 1-1 


OUTPUT results 


-9IT» U»mD» I 

UlO FOXMAT t//»99A»64HTi>4f following ST06>* SPECOS AND OJPECTIOXS GAVE A 
ISMEAX I\OE* CF ini//l 
w«|TE (6*m2i 

1U2 FOMVat I2IA»*ST0Py speed* • 2X»l9HSTUk'A OlPECTlONijA t J hDEV#6X» ** 

I lfHSJ^F-A90 A^ t9AGE*AXalSHS9D-)0C AvESAGf » 

9ITF I6«li:4) (^‘SS2SM*J>tSS2Sn*JUysU2S( Uji*S02SM*JltDEV2S 

I I l*Jl tS9wS2S( t *J» •S'^-'>2SM fJttSPO'»IOt I • »0I I 

111- »AT :2*x*r ' ^'SOk«:2*:a«*' 6«2*'«'* TEG» JA*A5t5X# 

IFS#2*4h v< S tlx 6*2 *4)^ DIG* ix«F b»2*4- ‘VSt^ AtFAt2*-<^ PEOI 
GO TO R20 

n A91TE (AtllUt I 

Ills »G«’i^AT (//tiixt • rHE‘<f wc^r so sto^m speeds a\"» dkectio^s that 

IGAvr A ShTAX ISGEA OF**lP 

•2r> c:st;nu' 

CD TO PJC 

92D aIITF UtlllRl 

;ilA FC^VM |//•A^•X•)6MT•4tS SnuNDISG NAD A SmCAE iSOfX OF 0*//l 
ttSITf lAtUWI 

»XITE tA»UUI ISAVEStSAVf StlSAvCOiSAvEDtSAVEUVtSAvfcXStSAVEwOt 

I S^ViWOtSAvE^S 

«S9 w 9:TC tAtlKr*) EUlSDSMPtESl 
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nZQ fO»MAT (//.47 x«SmE 1> •r6#2»7n# SI « •U«12 h* A‘<D €SI • •rs«2//l 

lOS CAlU :».::cX!TTIfSAt'AT#S30r«£ft»O«*Ml5S» 

S07 If (F'?'?CR.£0*'‘‘. .and. so T'J 1220 

' 10* IF !F'»ROR*t0.C. » iv‘R:»eC6fll541 TTUSwCAT 

309 U94 FORMAT 117X#*TTI « S^lO.l# •S*CAT • 

310 1202 IF (EORCR «F0. 1.1 W9ITE (6«11261 

- 311 1204 !► (ERROR .EO. 2*1 mRlTE (6fll2d) 

312 1206 IF (F9R0R .EC. S.l M91TF (Atll30) 

313 1126 FORvaT a7Xt*S^£AT TA'r.Cr 9E CALCuUTEO ^:ISSIN5 TE*>|F»/I 

• 314 - -U29 FORMAT (1?k#» 5»£AT CA'«^0T Sr CALCUtATEO MISSIES w$ V) ~ 

315 115" FORvaT (17x»»SwFAT CANNOT B£ CALCULATED FISS1.\5 fiO */) 

316 1^ (^'ISStEO.ClAHlTE (6fin6) S^CT 

317 1156 FO-iVAT (l4Xt • SPOT « '*FlO,ii . ... 

318 JF (vlss .£0.1 1 aRIT£cS*U32» 

319 1157 FORVA? (17X# ‘SPOT CANNOT 2E CALCULATED '"ISSInO OATA*/) 

- -370 - GO TO 1225- ... 

371 1270 •RITE I6*ll4r'» TTJ .S^EaT •SPOT 

322 1140 FCr^AT (l/.x#‘TTI • 'iFlO.lt *S**EAT • SFlC.l* »AN0 SPOT • ••FlO.l) 

373 GO TO 1.75 

324 1274 WRITE (6.6220) II 1 .<.L .LL^t T vO.Tm«S,SI NSS 

525 6220 FC7VAT (4l6#3a5*7l 

- 176 .-1225 CO TO 100 - — 

32T END 




SUPROJTINE ^^OEK(TT!♦S.'.EAT.SPOT.£aSCR^^:IS£l 
D|vf\SiCN T»3). 7Dt21. a’F13> 

COMWCN TEMPO UOTE'<APO) .40(41 ltWS(4l> .hour* NDA7. P 
CO-^'CN I ST A 
ISTA 

THIS PRCSRAY calculates TmK SOFT S*EAT AND SPOT INDICES 
TD(1) • SURFACE DEtPOOT F 

T’^j2; * pc:n7 c 

TOO! •lOCMP* 0£4 AOINT C- . ..... 

Till « SURFACE TEVPERATuRE F 
T(2» • S3C va TEVPERATuRE C 
TOI ■ 50CV3 TEMPERATURE C 
mm • surface rti o speed in knots 

#;I21 • 850 MP ^IND 5C££0 

WO I • 500VP WINN SPEED - . - - 

/VA • 4lSD vEfPOG ANGuE 850 TO 50C MB 

wO(l) • SURFACE .WND DIRECTION 

.vD(2J • "50 vn WIND 5H2CT20N .... 

WDOl •5'>C VP DIRECTION 

▼ t;« total totals roEx 

ASlN« alTIvftcr setting IN INCnES 
^fT:^►"C0=•(P♦•\^♦A/T(.')l 
C0^« I aoi3.75«*n.l902e)»0.:.vT65»/2£H.O 
1.0/J.19C2S<* 
fact* 29.P21/I013.25 
fRRLR *0. 
v|SS *0 

CONVERT TO P ROC'S R JMTS 

|» ( lEvPf:) .r.R. TfY0Ol,Eu*99.9) CJ To 165 

IF (DT‘'P(2i.c .JR. DtfYPOI.EJ 9.9) GO TO 165 

IF ( .J^.'.%5(r2i.ii:.;».9) go to itt* 

IF 1-.VT( ^1 59 i.y .(M. aO(‘J1.CO.V?.9I Gl TO 167 

IF ( wriC72J.F0.9‘M.c .OR. WD( ?7 » .•‘C.99.9 ? GO TO 167 
TUI* Tfvn(ii#i,H ♦17.0 
T(2)P TCvri7) 
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346 

34 -» 

34 * 

349 

330 

331 ~ 

332 

333 

334 


333 
336 - 
357 

338 

339 

340 
361 
36 r 

363 

364 
363 
166 
367 
36 n 
360 

370 

371 

372 

37 i 

373 

376 . 

377 

378 

379 

380 
3“1 


38? 

383 

3«4 

383 

186 

387 

388 

389 

390 

391 
302 

373 

394 

393 

396 

397 





TM»* TFM*(3I 
K(l)« 4 S( 1 ) 4 l «94254 
V.<J?) » fciSC BJ •1*94254 
• ^ 5 ! 22 l^l* 94?54 
TDll) « 32*0 

T3r?»» OTtw^PC?} 

T0I3>« 0TFMPI3) - — 

CALCULATE TTI 

TTI« T(2J ♦TD(2» -2.0»TC3I 

NO NEGATIVE TEP*^S 0CCU9 IN SWfAT ro9ruLA 

IP CTTI *lT. 49*DJ TT1«49*0 
!P <T0I2) •L'^. C*0) T0I2U 0*0 

CAUulATP SOFT S-EAT - * 

C |P either ‘4<<2) O^i a<C3I less TmAN 15 ♦ FVA«C 

1^ C «..<(? 1 .LT*!**;! *03* (w<l 3) .LT.15*OnP^A*J.O 
IP *L^* 130. 0» .09. (^018 UCT. 250. OM SC TO 150 

1^ MW0(?2UL7. 210*0) *0<3. (aOI2Z'*GT* 310*0))G0 TO 15C 


/.VA 

1 

» aO(22)- 

W0(8) 




IP 


hVA 

*LT* 

180*0) 

FVA 


*2606 

2 P 


aVA 

*LT. 

140*0) 

FVA 


•40 

IF 


a’VA 

*LT. 

120.0) 

FVA 


• 4463 

IF 


.sVA 

.LT. 

110*0) 

FVA 


• 5511 

IP 


aVA 

.LT. 

100*0) 

FVA 


.6394 

{r 


■•VA 

.LT. 

90*0) 

FVA 


*9404 

IF 


WVA 

• LT* 

8x7,0) 

FVA 

• 1,0 - 

IF 


.VVA 

.LT. 

70*0) 

Fv'A 

• 

*9404 

IF 


NVA 

.LT. 

60.0) 

PVA 

a 

.6894 

IF 


4VA 

• LT. 

50*0) 

FVA 

• 

• 4413 

JP 


tsVA 

.LT* 

40*0) 

FVA 

9 

• 2340 

IF 


V.'VA 

.LT. 

30.0) 

FVA 

m 

.0937 

IF 


X'V A 

.LT. 

20*0) 

FVA 

m 

• 0532 

IF 


•'VA 

• LT. 

10.0) 

FVA 

9 

• 0 


TO 1*0 






180 FvA«0.0 - - - • 

C CALCULATE SWEAT 

160 S VFAT» 12.0»T0I2) ♦ 20*04CTTI-n9.0l 42.C* (2 ) ♦iUAl 3 »♦ 125*0»PVA 
30 TO 93 . 

165 nso9«i, 

50 TO 95 

166 C 8 P 0 R * 2 * - . .... 

GO TO 95 

167 E990R • 3* 

C CAuCUtATr SPOT 

c calculate autt'eter setting 

95 IF I TE*'*P{ 1 ) ,r3.99*9 *0t<* OTE^P ( 1 » .t:o.99*9 ) uO TO 250 
2 F < p( i ) •E:*99.9 .on. Hf 1 ».fJ. 99.91 vO TO 250 
t {:).':.99./ .Of-* ,vDt D.E^.99.9) 00 TO 250 

CORR* I l^(COP»Ht 1 ) /Pi 1 1 ••*190234) )•*5^^ 

AS1N*JFACT *Pr)»CORR) 

ALTT ■ lCO.«t30*-ASlNI 

IF(ITIl) .LT* ‘O*.)) *ANO* ASIN .LT.29.5) ALTT* .SMLTT 

IP fwr^m *gt. ?60*o) wst* ?»w8ni 

!P U0(1) *lE* 260.0) GO TO lOl 
GO TO 200 

101 IF (.0(1) .LI. 230*3) CO TO 110 
IF (TOC) *LT* 55.0) wST«-?*a<(1 I 
IP (TOC) *lT* 60*0) wST --A^d) 
wsr«M<( 1 ) 

CO TO 200 

lie IF (wOC) ^Lt^ ?uu*0) GO TO 120 


jkIGINAL page ISI 
OF POOR QUAUTE 
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NL 


39^ 

399 

400 

401 - 

402 

403 

404 

405 

406 


407 
40 Q 
405 

410 

411 

412 
wl3 

414 

415 

416 
*•17 

-418 


IF iTOfl) .lT* 55.0) i<ST«0 

If (TD(i) .lT. 60.0) W$T«W<(;j/2«C 

i»ST«iif<(l) 

00 TO 200 - - — 

120 IF llr.'Om .LE. 140. 0) GO TO ISO 
nrST«W<tl) 

IF ( Toil) .LT, 55.) GO TO 200 — 

IF (NOAY .\r.74) GO TO 125 

IFC f <N0.t0.4?9).A O.!M0U‘^.EQ.515n.0fi.(tNC.10.3H9).A\D.I tMOJR.EO. 

11115). 0^. EO. 2315M ») GO TO 2CC ’ — 

125 nST ■ aXUI# 2. 

GO TO 200 

130 IF (v«(Oa) .LT. 70,0) GO TO 140 - - 

WST«W<I 1 J 
GO TO 200 

l40 »$T*0.0 — -- - - - 

200 SPOT* (Tll)-60.0) ♦ <T0(lJ-55.0) ♦ ALTT ♦ jtST 
GO TO 225 

2 SO y.i$s * I * - * - > 

225 CONTINUE 
RETURN 


sentry . 


no 






Table 3 


Variable 

Computer value 

Calculated value 

Storm speed 

5.45 

5.44 

Storm direction 

274.63 

274.64 

SURF-850 AV.REL. SPEED 

8.10 

8.105 

SURF-850 AV.REL. DIR. 

179.97 

182.03 

550-300 AV. REL. SPEED 

9.92 

9.90 

550-300 AV. REL. DIR. 

315.47 

315.51 

Shear Index 

5 

5 
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